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搞 要 : 油 藏 孔隙 分 布 具有 分 形 特征 ， 如 何 将 分 形 描述 方法 应 用 到 油 藏 多 孔 介质 内 流体 的 输 运 特性 研究 过 程 中 ， 
目前 的 研究 方法 存在 一 定局 限 ， 主 要 是 没有 考虑 量 纲 统一 原则 和 尺度 效应 等 。 本 文 基于 量 纲 统一 原则 ， 利 用 分 形 
者 述 方法 , 对 油 藏 多 孔 介质 内 输 运 过 程 渗透 率 和 渗流 流量 进行 了 推导 , 得 到 了 输 运 特性 参数 的 单位 尺度 转换 系数 ， 
据 此 可 以 实现 基于 微观 分 形 特 征 描述 来 分 析 宏 观 输 运 过 程 。 本 文 建立 的 分 形 输 运 描述 方程 能 够 实现 不 同 尺度 层级 
司 的 转换 ， 得 到 的 宏观 特性 参数 结果 更 接近 实测 值 。 基 于 本 文 所 建 模型 ， 利 用 某 油 藏 数据 进行 了 分 析 验 证 ， 得 到 
了 输 运 流量 及 渗透 率 变化 规律 与 实际 油 藏 的 分 布 规律 
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Study on the fractal transport characteristics of reservoir pore 
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Abstract: The reservoir pore distribution has fractal characteristics. How to apply the fractal description method 
to the transport properties of the fluid in the porous medium of the reservoir, the current research methods have 
some limitations, which is mainly without considering the principle of dimensional unification and scale effect etc. 
In this paper, based on the principle of dimensional unification, the permeability and seepage flow in the porous 
media are deduced by using the fractal description method. The unit scale conversion coefficients of transport 
characteristics are obtained, which can be used to analyze the macroscopic transport process basis of microscopic 
fractal characteristics. In this paper, the fractal transport equation can be used to achieve the conversion between 
different scale levels, and the results of the macroscopic characteristic parameters are closer to the measured 
values. Based on the model established in this paper, the data of a reservoir are analyzed and verified. The 
variation law of the transportation flow and the permeability is consistent with the distribution rule of the actual 
reservoir. 
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5 间接 地 反映 出 地 下 油 藏 的 情况 ， 据 此 可 以 理解 和 预 
0 51 号 测 地 下 油 藏 动态 , 进而 达到 提高 油气 采 收 率 目的 中 
所 谓 多 孔 介 质 的 输 运 是 指 在 孔 附中 的 液体 、 气 对 孔 隐 结 构 的 研究 可 追溯 至 20 世纪 20 年 代 ， 


体 、 电 、 热 、 声 、 波 等 在 多 孔 介 质 里 的 传递 由。 有 由 NuttingBI 首 次 通过 毛细 管 曲线 理论 结合 铸 体 薄片 
关 多 和 孔 介质 中 流体 输 运 特性 的 研究 通常 与 油 藏 开采 的 方法 研究 储 层 的 孔隙 特征 。 此 后 ， 人 们 逐渐 意识 
相关 联 ， 主 要 在 于 油 藏 储 层 介质 的 特性 可 以 直接 或 到 实际 油 藏 的 孔隙 结构 相当 复杂 ， 用 传统 的 数学 基 


础 方法 并 不 能 很 好 地 描述 出 真实 的 油 藏 孔隙 概况 。 
20 世纪 80 年代 B.B.Mandelbrotd 在 总 结 了 自然 界 中 
非 规则 图 形 后 ， 首 次 提出 分 形 这 个 概念 ， 可 以 用 来 
处 理 极 不 规则 的 复杂 形体 。 此 后 ， 一 些 研究 表明 储 
层 岩 石 的 宏观 和 微观 空间 都 具有 良好 的 分 形 特 征 。 
Avnir 等 吕 利 用 分 子 吸 附 法 证 明了 储 层 岩 石 的 孔 阶 
结构 具有 分 形 性 质 ; Katz 和 Thompson 等 63 利用 扫 
描 电 镜 观 察 岩 石 断面 ， 发 现 各 种 砂岩 、 页 岩 及 碳酸 
盐 岩 在 3 到 4 个 数量 级 的 范围 内 是 自 相 似 的 ， 可 用 
分 形 维 数 来 描述 孔 隐 的 结构 特征 , 其 值 在 2.55~2.87 
之 间 , 统计 自 相似 的 尺度 范围 在 0.2 上 Am 到 50 [1 之 
间 ， 之 后 他 们 通过 统计 分 形 的 方法 预测 了 岩石 的 孔 
隙 度 的 值 ， 证 明 与 实测 值 基本 吻合 ，HansenI9 通 过 
盒 维 数 得 出 砂岩 的 体积 分 形 维 数 和 表面 分 形 维 数 以 
及 分 维 与 岩石 渗透 率 的 关系 ， 证 明和 岩石 的 渗透 率 依 
赖 于 分 形 维 数 的 大 小 。Kim Tae HU9 等 基于 分 形 理 
论 开发 出 分 形 离散 裂缝 网 络 代码 ， 建 立 FDFN (分 
形 离散 网 络 ) 模型 ,获得 油 藏 孔 际 度 值 ,， Xie SY 等 
0 利用 电子 显微镜 呈现 二 维 图 像 ， 及 盒 维 数 和 多 重 
分 形 参数 计算 分 形 维 数 ， 建 立 它们 与 油 藏 孔隙 度 和 
渗透 率 的 关系 ， 评 估 不 同 地 层 的 油 藏 储量 。 

随 着 科学 技术 的 进步 和 油气 勘探 不 断 发 展 深 
入 ， 国 内 的 学 者 也 在 努力 探索 中 。 同 登科 等 02 使 用 
分 形 维 数 与 分 形 指数 来 描述 分 形 油 藏 的 渗流 特征 ， 
分 析 了 分 形 参数 对 分 形 油 藏 压 力 的 影响 ， 得 出 的 
预测 结果 与 实际 油 藏 的 原始 资料 很 接近 ; 郁 伯 铭 等 
由 采用 分 形 理论 和 方法 建立 油 藏 多 孔 介质 的 分 形 渗 
流 模 型 , 推导 各 参数 的 分 形 计 算 公式 , 得 出 解析 解 ; 
刘 俊 亮 等 (人 根据 分 形 基 本 概念 ， 利 用 Sierpinski 
就 分 形 多 孔 介 质 的 孔隙 率 、 渗 透 率 等 特征 参数 进行 
推导 ， 讨 论 了 不 同 分 形 模型 下 渗透 率 一 分 形 维 数 的 
关系 与 它们 的 差异 ; 孔 祥 言 等 人 (4 对 双重 分 维 的 多 
孔 介 质 导出 了 油 藏 渗流 速度 、 渗 透 率 及 孔隙 率 的 公 
式 ， 绘 制 了 油 藏 地 层 的 样板 曲线 用 以 分 析 油 井 ， 王 
自明 等 5 把 碳酸 盐 岩 储 层 裂缝 的 分 形 特 征 与 渗透 
率 数据 结合 ， 建 立 了 该 储 层 渗透 率 的 变异 函数 ， 并 
证 明 与 实际 情况 吻合 ; 刘波 09 等 利用 压 示 毛管 压力 
数据 得 出 低 渗 透 油田 储 层 岩 石 的 分 形 维 数 ， 并 与 其 
他 描述 岩石 非 均 质 性 的 参数 进行 对 比 ， 得 出 平均 分 
形 维 数 越 大 ， 多 孔 介质 的 非 均 质 性 越 强 ， 油 田产 量 
越 低 、 含 水 率 越 高 ， 于 本 福 等 人 0 建立 了 油 藏 压力 
分 布 的 分 形 渗流 模型 ， 推 导出 动态 储 层 压 力 的 预测 
公式 ， 得 出 分 形 维 数 对 储 层 压 力 的 影响 。 
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综 上 所 述 ， 基 于 分 形 理论 来 研究 油 藏 多 孔 介 质 
内 流体 的 输 运 特性 问题 ， 研 究 者 们 主要 是 通过 建立 
分 形 油 气 藏 渗流 模型 ,近似 描述 油 藏 地 质 孔 隙 结构 ， 
旦 取得 了 一 定 的 科研 成 果 。 但 是 实际 孔隙 结构 的 
内 部 属性 和 空间 分 布 都 存在 不 均匀 的 变化 ， 决 定 了 
多 孔 介质 的 结构 特性 在 空间 不 同 的 方向 上 会 存在 变 
异性 ， 如 一 些 大 的 孔隙 尺度 为 毫米 或 厘米 层级 ， 而 
微小 孔隙 只 有 纳米 或 微米 层级 那么 大 ， 因 此 孔隙 的 
尺寸 会 跨越 几 个 尺度 层级 , 从 分 子 尺度 到 微观 尺度 ， 
或 者 是 从 微观 尺度 到 宏观 尺度 上 的 差异 。 虽 然 储 层 
孔隙 在 不 同 尺度 范 围 内 都 具有 分 形 特征 ， 但 对 油 藏 
输 运 特性 的 影响 区 别 很 大 ， 即 大 尺度 下 表现 的 是 
藏 孔 际 整体 的 性 质 和 规律 ， 而 小 尺度 则 是 局 部 的 体 
现 。 这 使 得 简单 的 将 两 者 结合 来 分 析 油 藏 会 存在 分 
形 尺度 效应 问题 ， 即 在 一 个 尺度 层级 下 总 结 的 规律 
在 另 一 个 尺度 层级 并 不 一 定 适 用 ， 从 而 使 计算 结果 
与 实际 的 值 会 有 较 大 的 偏离 。 鉴 此 ， 本 文 基于 量 纲 
统一 原则 ， 利 用 分 形 描述 方法 ， 推 导 油 藏 多 孔 介 质 
内 输 运 过 程 模型 方程 ， 并 进行 实际 分 析 。 
1 分 形 渗流 一 般 描 述 方法 

多 孔 介质 的 输 运 特性 是 由 一 系列 特征 参数 来 
述 的 ， 其 中 最 重要 的 两 个 物性 参数 是 孔隙 度 和 渗透 
率 ， 因 此 对 这 两 个 量 的 分 析 是 研究 者 们 关注 的 重要 
问题 之 一 。 渗 透 率 、 孔 际 度 及 其 关系 式 决定 于 多 了 筷 
介质 的 微观 孔隙 结构 ， 经 过 人 们 研究 发 现 ， 天 然 形 
成 的 多 孔 介 质 和 部 分 人 工 制造 的 多 孔 材 料 的 微 结 构 
具有 分 形 特征 。 借 助 分 形 几 何 理论 ， 提 出 了 多 种 理 
论 模型 来 描述 多 孔 介 质 的 微 结 构 ， 建 立 多 孔 介 质 输 
运 问题 与 分 形 维 数 的 关系 。 

渗透 率 与 孔 阶 率 轨 、 分 形 维度 D, 等 参数 存在 
一 定 的 函数 关系 ， 关 系 式 表示 如 下 

K=K(g,D,.…) (1) 

根据 流体 的 流动 状态 ， 利 用 毛细 管束 模型 推导 

得 出 渗透 率 的 表达 式 为 5 


ND,{1_ 5 2? ) 5 
二 中 | ms (2) 
8r (6 -1/)3-D;y 


式 (2) 中 ，N 为 单元 体 多 孔 介 质 的 弯曲 毛 引 
管 或 孔隙 的 总 数 。 定 义 += 工 /LL 为 多 和 孔 介质 的 平均 
迁 曲 度 ， 即 毛细 管 或 孔道 的 实际 长 度 与 宏观 长 度 的 
比值 。5= i /hax 为 孔 介质 中 最 小 和 最 大 孔隙 半 
径 之 比 。 在 该 方程 中 , 对 总 数 Y 求 取 中 的 一 些 系数 ， 
没有 明确 概念 ， 也 无 从 知晓 怎样 得 到 ， 对 平均 迁 曲 


度 z 定 义 不 明 确 ， 此 推导 不 够 精确 。 

第 二 种 模型 是 多 孔 介质 渗流 力学 模型 ， 从 渗流 
力学 的 角度 出 发 ， 利 用 Fanning 摩擦 系数 了 与 雷诺 
数 R, 在 一 定 条 件 下 呈现 的 对 数 线性 关系 , 得 到 渗透 
率 表 达 式 [6 


ed a 
2c \2-Dr i™ -1 人 人 42 -1)3-D; 
G3) 

式 (3) 中 ，c 为 常数 ,可 以 由 摩擦 系数 和 雷 
诺 数 RR, 求 出 。 但 公式 需要 限制 条 件 ， 即 ,<5 的 情 
况 下 ， 一 方面 对 多 孔 介 质 中 渗流 速度 等 有 限定 ， 另 
一 方面 提出 的 相关 公式 并 未 与 实验 结果 进行 比较 ， 
无 法 得 到 验证 。 

对 渗透 率 的 研究 邦 伯 铭 等 中 通过 建立 分 形 分 析 
解 模型 ， 结 合 哈 根 一 泊 萧 叶 方程 ， 导 出 多 孔 介质 总 
流量 的 关系 表达 式 

MAmax 
0=- [a(WaN) 

Mi 


in 


和 AP Dy +Dr x| 1 | 和 om .0 
lag a BL DD as 

(4) 

式 中 ，AP 为 进出 口 两 端 压力 差 ， /4 为 流体 的 


黏度 系数 , 工 为 弯曲 毛细 管 或 孔 际 的 直观 长 度 ，D 
指 多 孔 介 质 的 分 形 维 数 ，D; 是 迁 曲 度 分 形 维 数 。 
由 流量 公式 结合 达 西 定律 得 到 渗透 率 的 函数 关 
系 表达 式 


Dr -1 D: 
人 上 和 5 
AP4 128 A 3+Dr 一 D/ 


式 中 ，A 是 某 个 单元 体 的 横 截 面积 。 

公式 〈5) 中 没有 经 验 常 数 ， 且 每 一 个 参数 都 有 
确定 的 物理 意义 ， 但 美中不足 的 是 推导 的 公式 中 都 
没有 考虑 单位 尺度 需要 统一 的 问题 。 即 不 了 解 孔 辽 
变化 规律 和 对 尺度 的 依赖 性 ， 只 是 单纯 的 自行 代入 
各 个 相关 参数 ， 随 意 取 值 ， 缺 少 理论 依据 ， 而 且 会 
产生 数据 量 纲 的 不 统一 问题 ， 这 势必 与 实际 产生 的 
结果 造成 偏离 。 因 为 渗透 率 、 流 量 等 特性 参数 是 描 
述 宏 观 多 孔 介质 的 ， 而 分 形 维 数 是 表征 微观 孔隙 结 
构 ， 在 已 知 某 一 尺度 下 的 性 质 ， 来 求 另 一 尺度 下 的 
性 质 时 ， 需 要 一 个 媒介 来 使 微观 尺度 的 描述 和 宏观 
尺度 计算 的 量 纲 达到 统一 也 就 是 本 文 研究 的 要 点 ， 
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在 考虑 单位 尺度 层级 统一 转换 的 基础 上 ， 对 以 上 提 
到 的 菜 一 公式 加 以 修正 ， 得 出 单位 尺度 转换 系数 ， 
减 小 误差 ， 以 更 加 接近 真实 值 。 


2 基于 量 纲 统一 原则 的 分 形 渗流 描述 
三 

前 述 描述 方法 存在 一 定 缺陷 。 首 先 ， 在 单 根 毛 
细 管 流量 的 表达 式 中 涉及 到 两 个 不 同 尺 度 层 级 的 参 
数 , 即 毛管 长 度 世 属于 宏观 尺度 , 毛细 管 管 径 人 1 是 
微观 尺度 ， 而 宏观 与 微观 对 孔隙 特性 的 影响 区 别 很 
大 。 其 次 ， 在 公式 推导 过 程 中 ， 并 没有 对 公式 中 每 
个 特征 参数 的 量 纲 尺度 下 定义 。 因 此 在 反映 孔 际 
中 流量 的 分 布 规律 时 ， 该 公式 的 右边 没有 保持 量 纲 
的 一 致 性 和 单位 的 一 致 性 。 根 据 这 一 原理 ， 此 渗流 
基本 公式 是 不 完整 不 完善 的 。 那 么 ， 二 者 结合 在 整 
个 公式 中 的 计算 时 ， 必 须 将 各 项 达到 量 纲 的 统一 ， 
需要 存在 一 个 某 尺 度 下 的 量 纲 转换 系数 。 

从 最 初 含有 长 度 工 和 管 径 人 的 方程 式 入 手 ， 
统一 为 宏观 尺度 ， 得 到 如 下 关系 式 


LY 
L -( 友 | KA=k PIP A ?=k PL (6) 


k 代 表 任 意 不 同 尺度 层级 间 需 要 达到 一 致 的 数 
量 级 ， 它 由 人 为 制定 的 尺度 与 国际 化 统一 量 纲 的 关 
系 确 定 。 

将 转换 后 的 毛管 长 度 方 程 代 入 流量 方程 可 得 到 


(= x AP (4 和 AP (x2) 
"0 Tap 0 hl 128 AZ 
AP (kA) 


128 HA kl ?ror A Yr 


=k3Dr A AP 入 
128 4 LA Yr 


单元 体 截面 的 流量 为 


=k*Pra(4) (7) 


Mnax 
0 =- a WaN) 

ii 
dp 和 AP 入 
全 128 7 六 和 新 


至 一 天 3+D7 


x DA0 pa 


和 3+D7 x AP 忆 / 


Mr =p PrO (8) 
128 1 IL? G+ Dr; -Di 


定义 K, = 2 ， 称 为 流量 尺度 转换 系数 。 
利用 Dancy 定律 ， 渗 透 率 的 公式 为 
ee (9) 
APA 
由 达 西 定律 可 以 知道 ， 公 式 中 的 相关 物理 量 都 
是 宏观 物理 量 ， 因 此 不 存在 尺度 层级 ， 直 接 代 入 流 
量 O ， 得 到 如 下 新 的 渗透 率 公式 
天 -人 三 于 Dy HLQ _ 3tprg (10) 
APA APA 
由 此 得 到 ， 分 形 渗透 率 公 式 的 转换 系数 
Kx=K,=k 。 
该 修正 公式 的 意义 在 于 ， 在 描述 宏观 物理 特性 
时 ， 优 先 考虑 某 量 纲 尺 度 下 数值 的 计算 ， 然 后 与 尺 


度 转换 系数 结合 ， 达 到 方程 两 边 量 纲 统一 的 目的 。 1 砂岩 Bl 的 三 维 扫描 图 像 
A 站 形 深 六 基 玫 太 we (像素 分 辨 率 5.34 LU ， 样 品 尺寸 2.138 mm) 
3 修 正 后 人 形 人 这 全 本 人 式 的 验 证 及 Fig.1 The three dimensional scanning image of sandstone B1 


应 用 


基于 本 文 建立 的 油 藏 多 孔 介质 的 流量 和 渗透 率 
模型 方程 ， 对 4 种 岩心 的 宏观 输 运 特性 参数 流量 和 
渗透 率 进行 数值 计算 ， 据 此 分 析 微 观 孔 际 特征 对 宏 
观 输 运 参 数 的 影响 规律 。 

3.1 实验 数据 及 基本 参数 

本 文采 用 4 个 不 同 油田 岩心 资料 进行 数值 计 
算 ， 并 将 修正 方程 的 计算 值 与 原 公式 以 及 实测 值 进 
行 比 较 。 这 些 样品 包括 3 种 砂岩 和 一 种 碳酸 盐 宕 ， 
其 原始 三 维 数据 是 通过 微 CT 扫描 成 像 系统 采集 得 
的 ， 能 直观 描述 岩心 特性 。 图 1 为 砂岩 B1 数字 
心 的 三 维 图 像 , 图 像 中 无 色 透 明 的 部 分 代表 和 孔 辽 ， 
色 的 代表 固体 骨架 ，X ,7,Z 分 别 代表 不 同 的 方 
; 图 2 为 B1 在 Z 方 向 上 二 维 截 图 所 得 到 的 图 像 ， ed 
其 中 白色 为 孔 险 ， 黑 色 则 是 固体 。 (像素 分 辨 率 5.34 jn ， 样 品 尺 寸 2.138 mm) 
通过 计算 岩心 的 基本 物理 参数 ， 从 而 更 加 精确 Fig.2 The two dimensional section of sandstone B1 
的 描述 分 析 油 藏 .实验 数据 及 基本 参数 如 表 1 所 示 。 

3.2 分 形 维 数 


(Pixel resolution 5.34 um , Sample size 2.138 mm ) 


志 


想 过 让 旦 
水 


(Pixel resolution 5.34 Wm , Sample Size 2.138 mm ) 


以 B1 岩心 为 例 ， 对 修正 公式 进行 验证 。 首 先 表 1 4 种 岩心 样品 实验 数据 及 基本 参数 
由 扫描 设备 得 到 所 需 数据 如 下 ， 最小、 最 大 孔 阶 直 Table 1 Experimental data and basic parameters of 
径 和 4, =5.345wm 、 ,=1256081m ， 孔 孙 度 4 Core Sarmples 
$=0.196， 渗 透 率 K=1286nD， 岩 样 宏观 长 度 由 样 号 二 SP 0 i 
=2.138gnm， 取 动力 粘度 =lmPars ， 压 差 mm {en EA YD 
AP=1Pa。 SandstoneB1 2.138 5.345 125.608 0.196 1286.0 
根据 微观 参数 ， 计 算出 孔隙 分 布 的 分 形 维 数 CarbonateCl 1.140 2.850 96.900 0.233 ”1102.0 


SandstoneS1 1.486 4.956 133.812 0.246 3898.0 
SandstoneS4 1.44 4.803 148.893 0.250 6966.0 


Dj 为 2.484; 迁 曲 度 分 形 维 数 D7 为 1.198。 男 外 ， 


其 余 3 组 岩心 与 此 相似 ， 不 再 累 述 。 
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透 率 的 值 ， 得 出 当 参数 取 值 在 不 同 尺度 层级 时 ， 计 
3. 3 修正 前 后 对 比 计算 分 析 算 结果 会 产生 不 同 程度 的 误差 ， 即 存在 尺度 效应 问 
根据 前 述 已 有 的 公式 和 修正 后 的 模型 ， 进 行 了 题 。 这 是 因为 油 藏 多 孔 介质 非 均 质 性 的 存在 ， 使 得 
计算 渗流 流量 和 渗透 率 的 对 比 计算 ， 详 见 表 2 及 表 空间 不 同位 置 上 有 不 同 尺度 大 小 的 孔 隐 结 构 ; 
3。 由 以 上 数值 计算 ， 可 以 得 到 以 下 结论 : 
(1) 通 过 比较 微米 尺度 与 毫米 尺度 下 流量 和 渗 


表 2 4 种 岩心 渗流 流量 对 比 
Table 2 Comparison of seepage flow in 4 cores 


流量 原 值 。” 微米 尺度 。 偏差 毫米 尺度 偏差 修正 流量 偏差 


lpm/s OQm 9% Qo py% 0 Gy% 
1 Bl 2. 484 2.36 14. 18 2. 33 15.27 2.765 -0. 55 
2 C1 1. 256 2. 04 -62. 4 2.037 -62. 14 1.752 一 39. 5 
3 Sl 5.792 5.677 1. 98 5. 667 2. 17 5.788 0.07 
4 S4 10. 03 8. 952 10.7 8. 944 10. 84 9. 524 5. 04 
表 3 4 种 岩心 渗透 率 对 比 
Table 3 Comparison of 4 kinds of core permeability 
册 渗透 率 ”微米 尺度 。 “偏差 ”毫米 尺度 偏差。 修正 渗透 率 偏差 
序号 上 
原 值 mD Own 111% Km 17127920 K 11731% 
1 Bl 1286 1105.1 14.07 1104 14. 15 1293 -0. 54 
2 Cl 1102 1789. 54 -62. 4 1782 -61.7 1536. 84 -39. 46 
3 Sy] 3898 3816. 42 2.1 3826 2. 00 3895. 02 0. 076 
4 S4 6966 6216.75 10. 75 6216.0 10.76 6613. 90 5. 05 
(2) 由 表 中 凶 号 岩 样 可 以 看 出 , 在 毫米 尺度 下 0 
的 值 比 微米 尺度 更 接近 实测 值 ， 但 是 得 出 的 结果 都 | 
与 实测 值 偏 离 很 远 。 分 析 原 因 ， 一 方面 可 能 是 提取 
的 岩心 样品 在 这 一 区 域内 孔 隐 尺寸 都 偏 大 ， 所 以 更 | 2000 
适合 尺度 层级 比较 大 的 运算 ， 另 一 方面 可 能 是 该 尺 让 全 
度 范围 内 的 此 种 分 形 特征 并 不 是 很 突出 ， 因 此 误差 二 1 
很 大 。 由 此 得 出 ， 虽 然 分 析 的 是 微观 孔隙 结构 ， 但 。 装 _ 1000 当 
是 ， 尺 度 范围 并 不 是 越 小 越 好 ， 而 且 在 研究 油 藏 多 
孔 介 质 时 ,单一 尺度 并 不 能 反映 出 油 藏 整体 的 特征 ; | 
(3) 对 比分 析 修正 前 后 的 渗透 率 及 流量 值 , 修 半生 生生 上 
正 后 的 结果 误差 更 小 ， 与 实际 值 更 接近 ， 所 以 修正 Wo 


公式 成 立 。 该 公式 的 修正 是 利用 方程 两 边 要 达到 量 
纲 统一 的 原则 ， 把 某 量 纲 微观 尺度 下 的 值 与 该 尺度 
下 的 量 纲 转换 系数 相 结 合 ， 通 过 系数 这 个 桥梁 ， 得 


3 孔隙 度 对 油 藏 流量 及 渗透 率 的 影响 


Fig.3 Effects ofporosity on reservoir flow and permeability 


出 最 络 宏观 结果 。 透 率 和 流量 的 变化 规律 及 影响 ， 结 果 如 图 3- 图 6 所 
3. 4 参数 变化 特性 影响 分 析 示 。 
选取 @ 号 岩 样 进行 参数 变化 特性 分 析 。 计 算 过 


程 中 通过 改变 岩 样 尺 寸 、 最 小 和 最 大 孔隙 直径 以 及 和 渗透 率 ， 随 者 孔 隐 度 增 六 而 增 六 。 这 是 因为 当 岩 
孔 阶 度 等 , 据 此 得 到 不 同 的 渗透 率 和 流量 的 变化 值 ，。 心 孔 辽 度 4 增 庆 ， 微 观 孔 隐 分 形 维 数 Dy 增 大 ， 由 此 
以 分 析 微 观 孔 隙 结构 的 参数 对 宏观 输 运 特性 参数 渗 得 到 孔隙 结构 非 均 质 性 增强 ， 和 孔隙 分 布 不 均匀 ， 会 
对 孔隙 内 的 渗流 性 能 在 不 同方 向 上 产生 一 定 的 影 


响 。 和 孔隙 度 $ 增 大 ,孔隙 平均 迁 曲 度 z, 减 小 ,， 迁 曲 


度 分 形 维 数 Dj 减 小 ， 筷 际 通 道 逐 渐 趋 于 一 条 直线 ， 
便于 孔隙 内 流体 流动 , 因而 渗透 率 和 流量 逐渐 增 大 。 
| F18Sx10? 
3.5x10° -| 入 
aouo*] LS Mis 
| ne 
2.5x10° -| Me 
仿 | A Te, FLA10 一 
和 2.0x10° -| ~ ee 在 
询 上 1zx10? 如 
于 和 
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oo 一 一 一 一 二 3.x107 
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图 4 最 小 孔隙 直径 对 流量 及 渗透 率 的 影响 
Fig.4 The influence of the minimum pore diameter on flow 
rate and permeability 

由 图 4 可知， 在 孔隙 度 一 定 的 情况 下， 渗透 率 
和 流量 随 着 最 小 孔隙 直径 的 增 大 而 减 小 。 同 理 ， 当 
最 小 孔隙 直径 4 逐渐 增 大 时 ， 微 观 孔 隙 分 形 维 数 
Dy 逐渐 减 小 ， 岩 心 内 不 同 大 小 尺度 的 孔隙 分 布 趋 于 
均匀 。 但 此 时 孔 际 的 迁 曲 度 分 形 维 数 Dj 逐渐 增 大 ， 
表征 孔隙 结构 内 流体 流通 通道 的 迁 回 曲折 程度 增 
大 ， 流 体 很 难 顺畅 的 穿 过 和 孔隙， 因而 此 时 流量 和 渗 
透 率 是 减 小 的 。 
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5 最 大 孔隙 直径 对 流量 及 渗透 率 的 影响 
Fig.5 The influence of the maximum pore diameter on 


flow rate and permeability 

日 图 5 可 知 ， 增 大 最 大 孔隙 直径 ， 岩 心 渗透 率 
和 流量 也 逐渐 增 大 。 因 为 在 孔隙 率 少 一 定 的 情况 下 ， 
最 小 孔隙 直径 Wan 不 变 ， 最 大 孔 阶 直径 4 不 

断 增 大 ， 由 此 计算 得 到 孔隙 分 形 维 数 Dy 增 大 ， 即 
岩心 非 均 质 性 增强 ， 不 同方 向 上 最 大 与 最 小 孔隙 之 
间 相 差 数量 级 逐渐 增 大 ， 此 时 微小 孔 隐 直径 可 以 忽 
略 ， 对 整体 流量 和 渗透 率 没 有 影响 。 并 且 随 着 0， 
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的 逐渐 增 大 , Dj 逐渐 减 小 , 孔 际 通道 趋 于 一 条 直线 ， 
便于 宏观 输 运 ， 因 此 渗透 率 和 流量 是 逐渐 增 大 的 。 
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抽样 尺寸 mm 
6 岩心 尺寸 对 流量 及 渗透 率 的 影响 

Fig.6 Effect of core size on flow rate and permeability 

由 图 6 可 知 ， 随 着 储 层 岩 心 取样 尺寸 的 增加 ， 
油 藏 的 流量 和 渗透 率 逐 渐 降 低 。 岩 心 直 观 长 度 尺寸 
石 的 增长 ， 分 形 维 数 Dy 并 没有 发 生变 化 ， 因 为 岩 
心 长 度 并 没有 对 岩心 的 内 部 结构 产生 影响 ， 因 此 此 
时 Dy 不 能 用 来 判断 其 对 宏观 输 运 特 性 参数 的 影响 。 
刀 增 大 孔 际 的 迁 曲 度 分 形 维 数 Dj 减 小 ， 但 因 岩 心 
两 端 压 差 不 变 ， 当 距离 增 大 ， 流 量 和 渗透 率 依然 减 
小 。 
4 结论 

多 和 孔 介 质 结构 组 成 极其 复杂 ， 精 确 描述 困难 ， 
为 此 一 些 研 究 者 们 尝试 引入 分 形 理论 ， 对 其 进行 
述 和 分 析 。 但 由 于 实际 多 孔 介 质 ， 和 孔隙 在 空间 不 同 
方向 上 存在 尺度 差异 ， 对 宏观 孔 际 输 运 特性 影响 很 
大 ， 尺 度 效 应 等 影响 不 应 忽略 。 本 文 考虑 量 纲 统 
原则 和 尺度 效应 ， 对 多 孔 介质 相关 分 形 输 运 模型 进 
行 了 修正 和 完善 ， 并 进行 了 计算 分 析 ， 得 到 以 下 结 


论 : 


(1) 针对 已 有 多 孔 介质 分 形 渗流 模型 存在 单 
位 、 尺 度 不 明确 问题 ， 本 文 进行 重新 推导 ， 得 到 了 
相关 的 分 形 输 运 分 析 模 型 ， 并 将 其 应 用 到 油 藏 多 孔 
介质 中 进行 了 验证 ; 

(2) 将 微观 孔 际 尺度 与 宏观 计算 分 析 相 结合 ， 
得 到 相关 输 运 特性 参数 ， 据 此 可 以 分 析 不 同 层级 间 
的 尺度 效应 ; 实际 分 析 表 明 ， 对 于 油 藏 并 不 是 尺度 
越 小 计算 结果 越 精确 ; 

(3) 基于 本 文 提出 的 分 形 输 运 模型 , 利用 某 阐 
田 油 藏 数据 进行 了 分 析 计 算 ， 并 分 析 了 微观 孔 际 特 
征 参数 对 油 藏 宏观 输 运 参 数 的 影响 ,实际 分 析 表 明 ， 
随 着 孔隙 度 和 最 大 孔隙 直径 的 增 大 ， 流 量 和 渗透 率 
也 增 大 ， 随 着 最 小 孔隙 直径 与 取样 尺寸 的 增 大 ， 流 


量 和 渗透 率 减 小 。 本 文 模拟 计算 分 析 结 果 与 油 藏 实 
际 输 运 特性 规律 相 吻 合 。 
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